The effect of spin-lattice relaxation on multiple-quantum NMR coherences in solids 
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A phenomenological theory of spin-lattice relaxation of multiple-quantum coherences in systems 
of two dipolar coupled spins at low temperatures is developed. Intensities of multiple-quantum NMR 
coherences depending on the spin-lattice relaxation time are obtained. It is shown that the theory 
is also applicable to finite spin chains when the approximation of nearest neighbour interaction is 
used. An application of this theory to an estimation of the influence of decoherence processes on 
quantum entanglement and its fluctuations is briefly discussed. 
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I. INTRODUCTION 

Multiple-quantum (MQ) coherences in a system of 
dipolar coupled spins can be detected by MQ NMR 
methods [1, 2]. In the simplest cases (one-dimensional 
spin chains and rings, a system of spin-carrying atoms 
(molecules) in nanopores), analytical and numerical 
methods were developed [1, Q in order to calculate the 
intensities of MQ NMR coherences. Usually spin-lattice 
relaxation is not taken into account at investigations of 
MQ NMR dynamics because the characteristic times of 
spin-lattice relaxation are longer than the characteris- 
tic time of MQ NMR dynamics [!]■ At the same time, 
MQ NMR methods have been also used to measure de- 
coherence for highly correlated spin states [1, Q and to 
study the scaling of the decoherence rate with the number 
of correlated spins @ . MQ NMR investigations showed 
both experimentally Q and with computer simulations 
for a system of equivalent spins Q that the relaxation 
time decreases with the increase in the MQ coherence 
order and in the number of spins. It was also shown Q 
that MQ NMR allows us to find a relationship between 
the intensity of the MQ coherence of the second order 
and a measure of the quantum entanglement which is 
the main resource of quantum computers [loj . However, 
that relationship was found without taking into account 
effects of spin-lattice relaxation [!|. A theoretical inter- 
pretation of such results is impossible without taking into 
account those effects. A theory of spin-lattice relaxation 
is also useful for traditional problems of NMR (paramag- 
netic impurities, internal molecular motions, etc.) jTl| . 

In the present paper we develop a phenomenological 
theory of spin-lattice relaxation for a system consisting 
of two dipolar coupled spins in the MQ NMR experi- 
ment at low temperatures. We calculate the intensities 
of MQ coherences of the zeroth and plus/minus second 
orders and discuss the effect of spin-lattice relaxation on 
MQ NMR dynamics. We show that the developed ap- 
proach is also applicable for linear spin chains when the 



approximation of nearest neighbor interactions is used. 
We discuss a connection of quantum entanglement with 
the intensity of MQ NMR coherence of the second order 
at different spin-lattice relaxation times. 



II. MQ NMR OF SPIN PAIRS WITH 
SPIN-LATTICE RELAXATION AT LOW 
TEMPERATURES 

We consider a system consisting of two dipolar coupled 
spins (s = 1/2) in the conditions of the MQ NMR exper- 
iment [1]. Initially, the system is in the thermodynamic 

equilibrium state in a strong external magnetic field i^o- 
The density matrix peq of this state in the laboratory 
reference frame is given by 



exp(ff4) 
Tr{exp(^4)} 



exp(^J.) 
4cosh^( 



(1) 



2kT^ 



where ojq = 7 is the Larmor frequency,7 the gyro- 
magnetic ratio, T the temperature, k the Bolzman con- 
stant, Iz — Izi + Iz2, and Iai{a. = x,y, z]i — \,2) \s the 
operator of the projection of spin i onto the axis a. 

We can neglect indirect spin-spin interactions and 
chemical shifts in Eq. ([Ij because they are much smaller 
than dipole-dipole interactions (DDI). The Hamiltonian 
of the DDI of the two spins, Hdz, in the strong external 
magnetic field fll| is 



Hd. = D{3hilz2 - hh), 



(2) 
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where D is the coupling constant and /1/2 = Ixilx2 + 

Iylly2 + Izllz2- 

We study MQ NMR dynamics of this two-spin system. 
The MQ NMR experiment [l| consists of four periods: 
preparation, evolution, mixing and detection. The spin 
system is irradiated on the preparation period by the 
multi-pulse sequence in which every period consists of 
eight resonance pulses [ij. It is suitable to switch to the 
rotating (with the Larmor frequency) reference frame [llj 
where the multi-pulse sequence consists of the four bloks 
of the form f/2-90° -2f-90^^-f/2 where f and 2f are 



2 



the time intervals between consecutive pulses and 90'^^ 
denotes pulses turning spins about axis ±x [l| . The DDI 
are oscillating fast, if the period 4*3-? = 12f of the multi- 
pulse seq uence is less than the inverse coupling constant 
As a result, spin dynamics of the system on 
the preparation period of the MQ NMR experiment is 
governed by the averaged non-secular Hamiltonian [l|, [l^ 



D 

Hmq ~ --^ 



(3) 



where and {i — 1,2) are raising and lowering op- 



erators (/j 



1,2). Spin dynamics on the 



preparation period is described by the Liouville equation 
for the density matrix p{t) (Tl| . 



* dt 



= [Hmq. pit)]. 



(4) 



and p{Q) — peq- It is worth to emphasize that MQ NMR 
coherences are determined by the density matrix at the 
end of the preparation period of the MQ NMR experi- 
ment di. 

In order to take into account spin-lattice relaxation 
one has to introduce, generally speaking, a non-linear 
term with respect to the density matrix p{t) in the right 
hand side of Eq. In practice, such equation can be 
solved with the linearization of the right hand side after 
it is expressed as an expansion in powers of p(t) — p, 
where p is the equilibrium density matrix in the RRF. 
This procedure is valid, if the deviation of p{t) from p 
is small. One can expect a small deviation, if the time 
independent p = p is a solution of the Liouville equation. 
It is reasonable to suppose that [l3{ 



exp(- 



hH 



MQ ' 



kT 



Tr {exp(- 



kT 



-)} 



(5) 



The temperature T is of an order of the lattice temper- 
ature. Its exact value is not significant in our considera- 
tion. 

We will consider the millikelvin temperature range, 
where one can suppose that |y ^ 1,^ = l/4e and 
e is the 4 by 4 unit matrix. Then the Liouville equation 
with a linearized phenomenological term for spin-lattice 
relaxation reads 



dp{t) 
dt 



-l[HMQ,p{t)] 



Pit) 



Tmq 



(6) 



where Tmq is the spin-lattice relaxation time of MQ 
NMR coherences. A reorienting spin group can be re- 
sponsible for relaxation 15]. We consider relaxation only 
as a pertubation of MQ NMR dynamics, i.e. we mean 
that D ^ 1/Tmq- We assume also that the spin-lattice 
relaxation time Ti in the laboratory reference frame is 
large (Ti Tmq) and we neglect it. It is worth to em- 
phasize that the equilibrium density matrix p describes 



the spin system in the RRF, where there is no any ex- 
ternal magnetic field, i.e. the Zeeman reservoir is ab- 
sent in the considered case. Thus, relaxation of MQ 
NMR coherences can be described with only one re- 
laxation time, Tmq- Introducing p*{t) — p{t) ~ p and 
p*[t) = p*(t)e*/'^Af« one can find that Eq. ^ can be 
written as 



d4^*{t) 
dt 



AHMQ,P*{t)], 



(7) 



and p*(0) = p(0) — p = Peq — P- At the end of the prepa- 
ration period r of the MQ NMR experiment, the solution 
of Eq. (O can be written as 



p*(r) = e-'^^^«>e9e'^^^«" - p. 



(8) 



The first term in the right hand side of Eq. (|8]) was ob- 
tained in 0, 



p-iHuQT iHmqt _ 

<= PeqC — ^ ^ 2 



4cosh^ 13/2 

(cosh /3-\-cos{Dt) sinh p —i sin(Z)T) sinh /3 \ 
8 J? S (9) 

2 siii(Z)r) sinh /3 cosh /3 — cos(£)r) sinh /3 / 

where 13= Using Eq. dH) one can find that 

p(r) = e-'""''^^ peqe'"'"<^^e-^/^'"<i + (1 - e^^/^*^«)p. 

(10) 

The density matrix pir + t) after the evolution period t of 
the MQ NMR experiment Q, during which the decoding 
field A (in frequency units) acts, is 

p{T + t) = ci-'^tl,^ e-'"MQr p^^^iHmqt ^tAtl, ^-t/Tmq ^ 

(1 -e"^/'^*^«)p. (11) 

We did not consider spin-lattice relaxation on the evo- 
lution period since we assumed that t <^ Tmq- Finally, 
after the mixing period t, during which the Hamiltonian 
Hmq changes its sign, the density matrix p(2r -I- 1) be- 
comes 

p(2t + t) = 

^iHMQr ^-lAtl^ ^-iHmqt p^^^iHmqt ^lAtl^ ^-iHmqt ^-2t /Tmq ^ 



(1 



-2r/TAfQ^- 



(12) 



Calculating the longitudinal magnetization 1^ (2t + t) at 
the end of the MQ NMR experiment [H and taking into 
account that Tr(/2p) = one can find that only MQ 
NMR coherences of zero and plus/minus second orders 
emerge in the our system and their intensities Jq{t), 
J±2{t), are: 

Jo(t) = tanh ^ cos2(i:»T)e~2^/'^^^« , 

J±2(t) = itanh^sin2(L>T)e"2r/T„Q_ ^^^g^ 
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Note that all performed calculations are also valid for 
any equilibrium density matrix p which commutes with 
operators Hmq and I^. 

It is also interesting to note that an analogous ap- 
proach can be applied to a linear chain of N interact- 
ing spins in the approximation of nearest neighbour in- 
teractions [T6l|. Only MQ NMR coherences of zero and 
plus/minus second orders emerge in such system. Using 
exact expressions for the intensities of these MQ NMR 
coherences , which were obtained without taking into 
account spin-lattice relaxation, one can obtain 



Jo(t) = 



-2t/Tmq 



y^cos^(2£'rcosfc), 



„ — 1t ITmq 

J±2(t) = ^sin2(2I?TC0sfc), (14) 



where k — {n = 1,2,..., N). It is worth to em- 

phasize that we neglect dipole-dipole interactions of next 
nearest neighbours which are only eight times weaker 
than the interactions of nearest neighbours. It means 
that the developed approach is justified only when 



D » 1/Tmq » D/8. 



(15) 



In considered systems we have found a simple exponential 
law for spin-lattice relaxation of MQ NMR coherences. 



III. QUANTUM ENTANGLED STATES IN MQ 
NMR 



It is well known jlO| that entanglement is the main re- 
source for quantum devices (in particular, quantum com- 
puters). MQ NMR is one of the numerous methods of 
obtaining entangled states [i,[l3. It turned out 9J that 
intensity of the MQ NMR coherence of the second order 
(without spin-lattice relaxation) is intimately connected 
with the concurrence [18,] which describes quantitatively 
entanglement of quantum systems. In order to calculate 
the concurrence jl8| at the end of the preparation period 
one should determine the matrix p(r) as 



/5(t) = {ffyl ® (Ty2)Pcc{T)i(Tyl (g) (T 



(16) 



where Pcc{t) is the complex conjugate matrix p(r) of 
Eq. (dni) in the standard basis {|00), |01), |10}, |11)} and 
ayi — 2Iyi, (i = 1, 2) is the Pauli matrix. 

The concurrence C of the two-spin system with the 
density matrix p{t) of Eq. pH)) equals to JJiJ 

C — max{0, 2A — Ai — A2 — A3 — A4}, 

A = maxj Ai, A2, A3, A4} (17) 

where Ai, A2, A3, and A4 are the square roots of the 
eigenvalues of the matrix product p{t)p{t). Although this 



product is not Hermitian, it has real nonnegative eigen- 
values [H. One can find from Eqs. (^1) . pU |) . ([TC ]) . ([T7 1) 
that the concurrence C(r) in our system is 



C(t) = max [ 0, 



''/Tmq 



2cosh^ /3/2 



I sin(£>T)| sinh /? — 1)- 

1 _ q-t/Tmq 



(18) 



This expression coincides with the one obtained in [qJ in 
the limit of the infinite spin-lattice relaxation time. The 
entangled state emerges at temperatures 



T <Te^ 



hujQ 



(19) 



where a = e^/^"«-l. If one takes wq = 27r500-10'^s"i the 
entangled state emerges at the temperature Te ~ 27mK 
when Tmq -^00^. According to Eq. the entangled 
states emerge at 2 > e^/^*^« - 1, i.e. at In 3 > t/Tmq- 
The temperature Te decreases with the decrease in the 
spin-lattice relaxation time. Using Eqs. (IT51) and ^TE\\ 
one can find a simple relationship between intensities of 
MQ NMR coherences of the second order and the con- 
currence: 



C 



tanh^[J2(r) + J_2(r)]- 



1 _ q-t/Tmq 



2cosh^ /3/2 



(20) 



This formula coincides again with the analogous one ob- 
tained in at Tmq 00. 

In order to obtain the entanglement measure of the 
two-qubit system we average over the state of the first 
spin and obtain the reduced density matrix for the second 
spin [1^. The reduced density matrix thus depends on 
the fluctuations of the first spin. These fluctuations giv e 
rise to the fluctuations of the entanglement measure [20| . 

The rms fluctuations, Ai?, of the entanglement entropy 
of the two-qubit system are [l^l 

Ai? = Clog2[^(l + v/l-C2)], (21) 



where the concurrence C is determined by Eq. ([Tl 

The dependencies of the sum of the MQ NMR coher- 
ences of the plus/minus second order, the concurrence 
C(t) and the entanglement fluctuations on the dimen- 
sionless spin-lattice relaxation rate t /Tmq at /3 = 6 are 
shown in Fig. [T] At long spin-lattice relaxation times the 
concurrence approximately coincides with the sum of the 
intensities of the MQ NMR coherences of the plus/minus 
second order as it was obtained in 0. At DTmq » 1, 
MQ NMR dynamics determines the entanglement and its 
fluctuations. Spin-lattice relaxation is the dominant fac- 
tor when DTmq ~ 1- We note also that the concurrence 
decreases with a decrease of the spin-lattice relaxation 
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FIG. 1: The sum of the MQ NMR coherences of the plus 
and minus second orders (sohd hne), concurrence C(t) (dash 
Une), the entanglement fluctuations (dash-dot line) versus the 
dimensionless spin-lattice relaxation rate t/Tmq- The spin- 
spin coupling constant D = 47rf307s~^, the duration of the 
preparation period r = 0.00086s, Dt = 97r/2, and /3 = 6. 



time and entanglement fluctuations are getting more im- 
portant. 

In conclusion we emphasize that quantum entangle- 
ment can be studied with MQ NMR methods and the 
effect of spin-lattice relaxation on MQ NMR coherences 
is important both for traditional magnetic resonance [llj 
and quantum information processing [lo| . 
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